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Relativistic electromagnetic penetration behavior is reexamined in the relativistic transparency region. The 
interaction is modeled by the relativistic hydrodynamic equations coupled with the full system of Maxwell 
equations, which are solved by a fully implicit energy-conserving numerical scheme. For the first time, we have 
studied the penetration behavior through ultra-short circularly polarized and linearly polarized laser pulses 
interaction with over-dense plasmas. It is shown that for the ultra-short circularly polarized laser penetration 
occurs through a soliton like behavior, which is quite consistent with the existing studies. However, we have 
found that the ultra-short linearly polarized laser penetrates through a breather like behavior, and there is 
an energy transition mechanism with a frequency 2co between the penetrated breather like structure and 
the background plasmas. A qualitative interpretation has been given to describe this mechanism of energy 
exchange. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw, 52.59.-f 



The dynamics of the penetration of ultra-intense laser 
irradiation into classically-forbidden plasmas has at- 
tracted wide attention for its distinct prospective appli- 
cations. The penetration of intense laser radiation into 
over-dense plasma core would be crucial in the concept 
of the fast ignition fusion 1 , plasma lens for ultra-intense 
laser focusing 2 and betatron oscillations for high intense 
electron acceleration 3 ^. Kaw et al. are the first to pre- 
dict that relativistically strong laser can reduce the ef- 
fective plasma frequency below the frequency of the laser 
radiation, allowing the laser to penetrate through the 
over-dense plasmas 5 -. This problem has been greatly en- 
riched by a series of theoretical and particle-in-cell simu- 
lation studies^—. Tushentsov et al. are the first to study 
the irradiation penetration of the circularly polarized 
(CP) laser pulse through fluid-Maxwell simulations 23 , 
but there is a serious shortcoming: the Maxwell equa- 
tions are solved within the framework of parabolic ap- 
proximation, which will certainly break down for ultra- 
short laser pulse. Afterward, Berezhiani et al. abandon 
the parabolic approximation and solve the full system 
of relativistic hydrodynamic and Maxwell equations 24 . 
However in their results, we have not found the intrinsic 
application to the study of ultra-short laser penetration 
into over-dense plasmas. In their fluid-Maxwell simula- 
tions, the laser duration and solitons formation time are 
over lOOXo- Their simulation results just reconfirm that 
of Tushentsov et al. with parabolic approximation. It 
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should be emphasized that their is no theoretical analy- 
sis on the penetration into over-dense plasmas with lin- 
early polarized (LP) laser, because the second-harmonic 
oscillating component of the ponderomotive force could 
disturb the background plasma density and make the 
problem even more complicated. The difficulty of the 
fluid-Maxwell simulations in studying LP solitons is that 
the explicit numerical scheme will break down in dealing 
with density ripples caused by second-harmonic oscillat- 
ing ponderomotive forced 5 .. However, the crucial issue, 
which is of fundamental influence to the fluid-Maxwell 
simulations, is the vacuum heating and J x B heat- 
ing mechanisms on the vacuum-plasma interface. These 
heatings are totally kinetic mechanisms which can not 
be described by a fluid treatment. Whereas in our work, 
these strait problems are solved in two ways. A fully im- 
plicit energy-conserving numerical scheme is applied to 
well deal with the density ripples caused by the oscillat- 
ing term. Abandoning the parabolic approximatio n 24 ' 25 , 
we solve the full system of relativistic hydrodynamic and 
Maxwell equations, and focus on ultra-short laser irradi- 
ation with full-width-half- maximum 5 To. In such cases 
the vacuum heating and J x B heating can be neglected 
considering the ultra-short interaction time between the 
laser and plasmas on vacuum-plasma interface. 

In this paper, the interaction between the ultra-short 
laser irradiation and the plasmas is modeled by the rel- 
ativistic hydrodynamic equations coupled with the full 
system of Maxwell equations, which are solved by a fully 
implicit energy-conserving numerical scheme. The pene- 
tration of the ultra-short CP irradiation occurs through 
a soliton like behavior which is consistent with previous 
theoretical and fluid-Maxwell simulation studies 23,24 ' 29 , 
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and its propagation in the homogeneous background plas- 
mas is quasi-static. However, it is shown that, for the 
ultra-short CP laser, the penetration occurs through a 
breaker like behavior, and there is an energy transition 
mechanism with a frequency 2o>o between the penetrated 
breather like structure and the background plasmas. We 
have proposed a qualitative interpretation to describe 
this mechanism of energy exchange. 

The governing equations are the relativistic hydrody- 
namic fluid equations coupled with the full system of 
Maxwell equations. As we are considering the relativis- 
tic intense laser irradiation and the electron quiver ve- 
locity under the laser electric field is much greater than 
the thermal velocity, we treat the plasmas as cold liq- 
uid. Below are the group of normalized full system of 
equation o 24 ' 29 : 
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7e = [l + a 2 + a 2 +p 2 ] 1 / 2 , (8) 

li = [l + {a y /mf + {a z /mf+p 2 \ 1 l 2 . (9) 

The group of equations are normalized in terms of dimen- 
sionless quantities t = coot, x = ujqx/c, a y = eA y /m e c 2 , 
a z = eAy/nieC 2 , n e = n e /n c , n, = n,/n c , p e = p e /m e c, 
Pi = Pi/m-iC, m — rrii/m e and <f> = ecj>/m e c 2 , where uiq is 
the laser frequency, a y and a z are the normalized ampli- 
tude of the laser irradiation, n e and rii are the normalized 
plasma density for electrons and ions, p e and pi are the 
normalized momentum along the laser propagation direc- 
tion for electrons and ions, m is the normalized ion mass, 
4> is the normalized electrostatic potential and n c is the 
critical density. 

To solve the group of equations, a fully implicit energy- 
conversing numerical scheme has been applied 2 ^, which 
can well deal with the density ripples and restrain over- 
shooting of the numerical scheme. Here we have assumed 
that the boundary conditions [a y — 0, a z — and = 
on both sides of the simulation box] are always satisfied 
during the simulation. At t — O.OTo the laser pulse is 
apart from plasmas, and at later time, the laser irradi- 
ates on the plasmas. Our simulation is a global fluid- 
Maxwell simulation, which is quite different from pre- 
vious studies that the soliton structure is put in as an 
initial conditio n 15 i 16 i 27 i 28 . The exact definition of a fluid 



implies the accumulation of a macroscopic number of par- 
ticles within the interest region, and it seems that the 
fluid treatment of the region of zero plasma density may 
result in unphysical outcome. When numerically inte- 
grating the equation of motion Eqs. (5) and (6) within 
the vacuum-like region, the lack of the balance force from 
the electrostatic field may result in the unphysical x di- 
rection speeding of the electrons and ions under the laser 
ponderomotive force. However, the plasma density is 
zero, and the total x direction current is zero too, which 
does not contribute to the solving of the coupled Maxwell 
equations and makes the global fluid-Maxwell simulation 
applicable^ 4 -. 

To ensure the accuracy of the numerical scheme, we 
follow the energy of the system at each time step during 
our simulations. The total energy can be separated into 
four components: 

,4/^) 2 + (^) 2 + (^) 2 + (^) 2 Kio) 

(11) 



dx[(-f e - l)n e 
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where e; is the laser energy, e p is the energy of the plasma 
wave, e e is the electron kinetic energy and is the ion 
kinetic energy. In Eq. (10), we can see that the laser en- 
ergy consists of four parts: y component electric energy 
{(da,y/dt) 2 ], y component magnetic energy [(da z /dx) 2 ], 
z component electric energy [(da z /dt) 2 ] and z compo- 
nent magnetic energy [(da y /dx) 2 ]. It should be empha- 
sized that this field is plagued by substandard simulation 
methodology, which has led to bold conclusions being 
drawn from what amounts to numerical artifacts. The 
energy conserving of the simulations convince us that the 
new physics we report is not from numerical effect. 

We first study the problem of an ultra-short CP laser 
pulse penetration into over-dense plasmas. The parame- 
ters of the incident CP laser pulse are: normalized am- 
plitude a y = 0.71, a z = 0.71, [sin(7r?;/2r) 2 ] time profile 
with r = 5.0To and laser frequency ljq = 0.88. The 
background plasma density is set to be n e — 1.0, cor- 
responding to n e = 1.3n c , in which n c = m e u) 2 /4ire 2 . 
The vacuum-plasma is modeled with [sin(7ra;/2s) 2 ] den- 
sity profile, where s is the shape factor to describe the 
steepness of the vacuum-plasma interface. 

In Fig. 1, the space-time dynamics of the penetration 
process is demonstrated when there is a sharp density 
slope at the plasma boundary. A relativistic electromag- 
netic soliton is generated during the penetration process, 
which is in agreement with the theoretical prediction that 
the penetration occurs through soliton-like structure 29 . 
Compared with the incident pulse, the penetrated soli- 
ton owns a sharply narrow width which is less than one 
laser wavelength and the amplitude is greatly enhanced. 
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FIG. 1. (color online) The formation process of CP solitons: red line represents the normalized laser intensity and the black 
line is for the plasma density. Here, a y = 0.71, a z = 0.71, r = 5Tb, cjq = 0.88, n e = 1.0 and s = 0.4Ao- 
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FIG. 2. (color online) The space-time distribution of the nor- 
malized laser intensity: (a) is for CP laser and (b) for LP 
laser. Here the simulation parameters for CP laser pulse are 
the same as that of Fig. 1, and the simulation parameters for 
LP laser pulse are a y — 1.00, a z = 0.0, r = 5Tb, ojq = 0.88, 
n e = 1.0 and s = O.lAo- 



After the formation, the soliton propagates into the uni- 
form plasmas with a constant velocity v = 0.172, which 
is shown in Fig. 2 (a). From Fig. 3, the energy history 
demonstrates that during the soliton formation stage, 
within the early 7.0Tb, there is an energy transfer from 
the incident laser to the background plasmas. It is easy 
to understand this energy transfer mechanism: the pon- 
deromotive force drives the plasmas inward in the laser 
propagation direction, during this process work has been 
done by the laser pulse. However during the propagation 
process of the soliton, we see from Fig. 3 that the energy 
of the laser pulse and background plasmas do not change 
any more, which demonstrates that the CP soliton prop- 
agation is quasi-static. 

Considerable part of the incident CP laser pulse is 
transformed into the sharply narrow and relativistic in- 
tense electromagnetic soliton. Our simulations demon- 
strate that there is an obvious relation between the pen- 
etration ratio and the steepness of the vacuum-plasma 
interface. As is shown in Fig. 3 (a)-(d), for a relatively 
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FIG. 3. (color online) The energy history of the fluid-Maxwell 
simulations: back line represents the normalized energy of the 
CP laser, the red line represents the plasma energy (mainly 
the electron kinetic energy) driven by CP laser, the green line 
is for energy of the LP laser and blue line is for the plasma 
energy by LP laser. Here the simulation parameters are the 
same as that of Fig. 2. 



steep interface the penetrated CP soliton is quite weak, 
and the penetration ratio is rather low. In contrast, the 
rather smooth interface allows a strong penetration. Fig. 
3 (e) shows the penetration ratio for different shape fac- 
tor s, demonstrating the general penetration ratio range 
10% -60%. 

Fig. 5 shows the spectral intensity of the initially in- 
cident CP irradiation, the reflected part and the pene- 
trated soliton. In the soliton formation stage, the pon- 
deromotive force pushes the plasmas inward in the direc- 
tion of the incident pulse propagation. As red line shows, 
the reflected radiation from the vacuum-plasma interface 
is mostly red-shifted. The frequency of the penetrated 
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FIG. 5. (color online) The spectral intensity of the incident 
CP laser pulse, reflected part and penetrated part. The black 
line is for initially incident CP laser pulse, red line is for re- 
flected and green line for penetrated, (a) is the case s = 0.1 
and (b) is the case s = 0.4. Here the remaining simulation 
parameters are the same as that of Fig. 1. 



FIG. 4. (color online) (a)-(d) demonstrate the reflected and 
penetrated CP laser pulse structure under different vacuum- 
plasma steepness, where the red line is the laser intensity and 
black line is the plasma electron density, (e) The penetration 
ratio of the incident laser pulse vs. shape factor s, where the 
black triangle is from the fluid-Maxwell simulation data. 



CP soliton is blue-shifted for s = 0.1 [subplot (a) in Fig. 
4] and red-shifted for s = 0.4 [subplot (c) in Fig. 4]. The 
behavior of the weak penetrated soliton [s = 0.1, subplot 
(a) in Fig. 4] can be well described by the known finite 
soliton velocity solutions, which are based on the weak 
density response^ 



(1 - v 2 - if) 



2\l/2 



rjv 



cl v = a m sech[ - (x — vt)] x cos[— '- — tt(x 

where 

. _ 1[ (l-v 2 )(l-v 2 -V 2 ) n /2 



4 V 2 _ (l-v 2 )(3 + v 2 )' 



is the amplitude of the CP soliton, and rj is the frequency 
parameter with i] = ui(l — v 2 ). The carrier frequency 
of the weak penetrated soliton can be derived from the 
frequency-amplitude- velocity relation, 



[16(1 - v 2 ) 2 + 4,(1 - v 2 )(3 + v 2 )} 1 / 2 
(l-w 2 )[4a^ + 16(l-w 2 )] 1 /2 



(14) 



For the weak penetrated soliton [s = 0.1, subplot (a) in 
Fig. 4] with a m = 0.482 and v = 0.172, Eq. (14) gives 
uj = 1.008 = 1.145woi which is quite consistent with the 
one 1.142wo obtained from the fluid-Maxwell simulation 
in Fig. 5 (a). 

To study the penetration process of the LP laser, we 
should be very careful. The density ripples caused by 
the second-harmonic oscillating ponderomotive force will 
break down the simulation process, while an fully im- 
plicit energy-conserving numerical scheme is applied to 
well deal with this problem. In addition, in order to 
avoid or neglect the vacuum-heating and the J x B effect 
on the vacuum-plasma interface, the incident laser pulse 
should be ultra-short to shorten the interaction time be- 
tween the incident laser and the sharp vacuum-plasma 



boundary. We change the incident laser amplitude (the 
same intensity with that of CP) to a y = 1.0 and a z = 0, 
while keeping other parameters the same as that of CP 
laser. Similar to the CP case, during the early 7.0Tq 
interaction time, the energy is transfered from the inci- 
dent LP laser pulse into background plasmas in one way 
although there is some small scale perturbations. Fig. 2 
(b) shows that the propagation velocity of the penetrated 
electromagnetic structure is almost the same as that of 
CP soliton which is v — 0.172. However, we find that un- 
like the CP soliton, during its propagation, there exist an 
energy transition mechanism between the LP penetrated 
structure and the background plasmas with a frequency 
about double the incident laser frequency. The energy 
)] i exchange is even more clearly indicated from the energy 
history map shown in Fig. 3, where the laser energy and 
the background plasma energy coupled together and ex- 
change with each other with a frequency 2ojq. 

However, from Fig. 6 (c) and (d) for the LP penetra- 
tion case, we can see that the penetrated electromagnetic 
structure disappears at t — 16.00To and suddenly reap- 
pears at t = 16.25T , with a period T /2. At t = 16.00T 
almost all the penetrated electromagnetic energy is trans- 
formed into background plasmas. While at t = 16.25To 
the electromagnetic structure has returned back from the 
background plasmas. It behaves like a blinking breather, 
which is quite different from the quasi-stable soliton as 
shown in Figs. 6 (a) and (b). 

Here, we provide a qualitative interpretation for the en- 
ergy transition mechanism with 2ujq between the LP pen- 
etrated structure and the background plasmas. During 
our simulations, we follow the kinetic energy of the elec- 
trons, the kinetic energy of the ions, the laser energy and 
the plasma wave energy. As the electron kinetic energy 
[e e ] occupies the main part of the background plasma en- 
ergy [ei + e e + e p ] , the plasma wave energy does not dom- 
inate during these process, thus we can drop the Poisson 
Equation [Eq. (7)] during our analysis. We now focus on 
the Eq. (1). The refraction term from the background ion 
motions can be neglected due to the high ion to electron 
mass ratio and we set the electron momentum along the 



5 



3.0 
2.4 
1.8 
1.8 
0.6 



0.0 

3.0 



3.4 
1.8 
1.3 
0.6 



fa; - 


w - 


: 1 : 


: A : 




/ A ■ 




\ A - 








A A 







-11 3 -SI 1 3 5 

xAo x/\ a 

FIG. 6. (color online) The distribution of laser intensity 
and plasma density: (a) and (b) is for CP laser pulse at 
t = 16.00T and t = 16.25Tb; (c) and (d) for LP laser pulse 
at t — 16.00Tb and t — 16.25To. The black line represents 
the plasma electron density, red line is for laser intensity and 
green line for electron kinetic energy density. Here the simu- 
lation parameters are the same as that of Fig. 2. 



laser propagation direction to p e 
be rewritten as 



d 2 a y d 2 a y 



0, then Eq. (1) can 



dt 2 dx 2 (l + ag) 1 / 



2 U J 



= 0. 



(15) 



As a y — a m sm(uiot), the refraction term n e /(l + 
a 2 ) 1 / 2 which describes the effective plasma density un- 
der the LP laser irradiation becomes n e /[l + — 
cos(2wo0] 1 ' /2 - I R our fluid-Maxwell simulations, 
the laser frequency is ujq — 0.88 which corresponds to 
n c = 0.76 and the background plasma density is n e = 1.0 
which corresponds to n e = 1.3n c . Unlike the soliton 
structure in the CP laser case where there is a density 
cavity for the background plasmas, in the LP laser case, 
the background density is almost unchanged, n e ~ 1.0, 
except for some small scale and fast oscillating density 
ripples. Taking into the simulation result, a m ^2.0, the 
effective plasma density n e /[l+a^/2— a^/2 cos(2uj t)] 1 ^ 2 
can range from 0.45 to 1.00 where the critical density 
n c located at 0.76. Let us consider the two extreme 
conditions. If the effective plasma density reaches 0.45, 
it means 7 e reaches maximum or the background plas- 
mas obtain its highest energy. However the low effec- 
tive plasma density can not capture the relatively high 
frequency penetrated structure with ujq = 0.88, and this 
penetrated electromagnetic structure just disperses away, 
as shown in Fig. 6 (c). If the effective plasma density 
arrives at 1.00, the background plasmas reach its low- 
est energy. The relatively low frequency electromagnetic 
structure is well trapped by the high effective plasma den- 
sity, which is shown in Fig. 6 (d). This process repeats 
at every half laser period as shown in Figs. 2 and 3. 



When the background plasma density n e is greater 
than 1.5n c , the interaction results in the generation of 
plasma field structure consisting of alternating plasma 
and vacuum regions with the CP laser penetrating into 
the over-dense plasmas over a finite length, which is con- 
sistent with the previous studiesSii^. Because the forma- 
tion time of the alternating plasma field structure is quite 
long which is about hundreds of laser period, it call for 
quite long incident laser pulse. As analyzed, it beyonds 
the ability of fluid-Maxwell simulations in studying such 
long pulse LP laser interaction with over-dense plasmas. 

In summary, the relativistic electromagnetic penetra- 
tion behavior is reexamined in the relativistic trans- 
parency region, which is modeled by the relativistic hy- 
drodynamic equations coupled with the full system of 
Maxwell equations. We have studied the penetration be- 
havior through ultra-short CP and LP laser pulses inter- 
action with over-dense plasmas. It is shown that for the 
ultra-short CP laser penetration occurs through a soliton 
like behavior, which is quite consistent with the existing 
studies. The penetration ratio of the ultra-short CP irra- 
diation is systematically studied with different steepness 
of the vacuum-plasma interface, which indicates that a 
rather smooth vacuum-plasma interface allows a strong 
penetration. The propagation of the CP soliton in ho- 
mogeneous background plasmas is quasi-static. However, 
we have found that the ultra-short LP laser penetrates 
through a breather like behavior, and there is an energy 
transition mechanism with a frequency 2wo between the 
penetrated breather like structure and the background 
plasmas. A qualitative interpretation has been given to 
describe this mechanism of energy exchange. 
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